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Forkhead box O (FOXO) transcription factors are key players in diverse cellular processes affecting 
tumorigenesis, stem cell maintenance and lifespan. To gain insight into the mechanisms of FOXO- 
regulated target gene expression, we studied genome-wide effects of FOX03 activation. Profiling 
RNA polymerase II changes shows that FOX03 regulates gene expression through transcription 
initiation. Correlative analysis of FOX03 and RNA polymerase II ChlP-seq profiles demonstrates 
FOX03 to act as a transcriptional activator. Furthermore, this analysis reveals a significant part of 
FOX03 gene regulation proceeds through enhancer regions. FOX03 binds to pre-existing enhancers 
and further activates these enhancers as shown by changes in histone acetylation and RNA 
polymerase II recruitment. In addition, FOX03 -mediated enhancer activation correlates with 
regulation of adjacent genes and pre-existence of chromatin loops between FOX03 bound enhancers 
and target genes. Combined, our data elucidate how FOXOs regulate gene transcription and provide 
insight into mechanisms by which FOXOs can induce different gene expression programs depending 
on chromatin architecture. 
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Introduction 

The Forkhead box O (FOXO) family of transcription factors 
consists of FOXOl, FOX03, FOX04 and FOX06. In model 
organisms, FOXOs function as tumour suppressors (reviewed 
in Dansen and Burgering, 2008), and affect both lifespan 
(reviewed in Partridge and Bruning, 2008; Kenyon, 2010) and 
stem cell maintenance (Miyamoto et al, 2007; Paik et al, 2007; 
Tothova et al, 2007; Miyamoto et al, 2008; Yalcin et al, 2008; 
Paik et al, 2009; Renault et al, 2009; Zhang et al, 2011). An 
increasing number of FOXO target genes has been identified, 
responsible for modulating a variety of cellular functions 
including cell death, proliferation, metabolism and scavenging 
of reactive oxygen species (reviewed in van der Vos and Coffer, 
2011). In general, two conserved pathways regulate FOXO 
activity. On the one hand, FOXOs are negatively regulated in 
response to insulin and growth factor signalling (Biggs et al, 
1999; Brunet et al, 1999; Kops et al, 1999; Nakae et al, 1999), 
while positive regulation occurs in the presence of oxidative 
stress through JNK activity (Brunet et al, 2004; Essers et al, 



2004; van den Berg and Burgering, 2011). Gene expression 
analysis of different cell types in Foxol, 3, 4 triple knockout 
mice has shown that FOXO target gene activation is highly 
context dependent. Little overlap in FOXO-regulated genes 
between different primary cell types was found, suggesting 
that FOXOs regulate target gene expression in a cell type- 
specific manner (Paik et al, 2007). This context-dependent 
regulation has been attributed to both the output of specific 
signalling pathways, resulting in different combinations of 
post-translational modifications on FOXO (Calnan and Brunet, 
2008) and to differential co-factor binding (van der Vos and 
Coffer, 2008, 2011). 

Combining genome-wide FOXO-binding profiles with 
analysis of FOXO-induced changes in gene expression is 
crucial for understanding FOXO biology. Studies in Drosophila 
and Caenorhabditis elegans (C. elegans) linked genome-wide 
binding of FOXO orthologues to changes in gene expression in 
FOXO-deficient strains and have revealed a complex process in 
which binding was not the sole determinant of gene regulation 
(Teleman et al, 2008; Schuster et al, 2010; Alic et al, 2011). 
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These studies are however limited, since changes in transcrip- 
tion are based on steady-state levels in mutants and likely 
include secondary and adaptive responses. Genome-wide 
binding studies of mammalian FOXOs using deep sequencing 
have been performed (Fan et al, 2010; Lin et al, 2010; Litvak 
et al, 2012), but are limited in the detection of direct 
consequences of FOXO activity. 

In this study, we use a mammalian cell culture system with 
inducible FOX03 to study direct consequences of FOX03 
activation. Changes in mRNA levels can be influenced by 
mRNA stability and might therefore not completely reflect 
direct effects of FOX03 induction. We therefore determined 
changes in RNA polymerase II (RNAPII) occupancy as a more 
direct measure of transcription factor activity (Mokry et al, 
2012). Combined with a FOX03 genome-wide binding profile, 
we show FOX03 acts as a transcriptional activator, regulating 
target gene expression through transcription initiation. We 
reveal a function for FOX03 at distal locations outside of 
promoter regions and show the existence of loops between 
several FOX03 bound distal regions and regulated genes. 
These data, together with the cell type specificity observed in 
other studies, suggest a model in which spatial organization of 
the genome is relevant for FOXO target gene regulation. 

Results 

Identification of FOX03-induced transcriptional 
changes through genome-wide profiling of 
RNAPII occupancy 

To study FOX03-mediated changes in gene expression, we 
used DLD1 colon carcinoma cells containing a 40H Tamoxifen 
(40HT) inducible FOX03A3-ER fusion (DLD1-F3 cells) 
(also called DL23 cells; Kops et al, 2002b). 

A constitutively active FOX03A3 mutant, lacking three 
inhibitory PKB phosphorylation sites, is fused to the ligand 
responsive domain of the oestrogen receptor. This domain 
contains a point mutation and is therefore transcriptionally 
inactive and only responsive to 40HT (Littlewood et al, 1995) . 
DLD1-F3 cells contain comparable levels of inducible 
FOX03A3 relative to endogenous FOX03 and induction with 
40HT results in the characteristic FOXO effects on cell cycle 
(Kops et al, 2002b) , making it a suitable system to study the 
global effects induced by the activation of a single transcrip- 
tion factor. Compared with physiological stimuli that activate 
endogenous FOXOs but also influence parallel pathways, this 
system allows us to study global effects induced by isolated 
activation of FOX03. We first determined RNAPII occupancy 
by chromatin immunoprecipitation combined with deep- 
sequencing (ChlP-seq) in DLD1 and DLD1-F3 cells, providing 
a snapshot of global RNAPII distribution before or after 4 and 
24 h upon addition of 40HT. We calculated changes in RNAPII 
occupancy at the coding region of all annotated genes. We 
observed increased RNAPII occupancy specifically in DLD1-F3 
cells after 40HT treatment in numerous genes and identified 
many known FOX03 targets as well as potential new 
FOX03 target genes (Figure 1A and B; Supplementary 
Table 1; Supplementary Figure 1). Cross-reference with genes 
regulated in micro-array analysis (data set recently published 
in Ferber et al, 2012) revealed considerable overlap, with a 
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higher degree of overlap for up- compared with downregulated 
genes (Figure 1 C) . This results in a list of target genes for which 
expression is altered through changes in transcription. 
Compared with FOXO micro-array studies, this excludes 
annotation of FOX03 targets through transcription-indepen- 
dent effects. 

To gain insight into the functional outcome of direct FOX03- 
induced changes, we performed a Gene Ontology (GO) term 
analysis. The subset of upregulated targets reveals enrichment 
of genes involved in signalling (Supplementary Table 1). 
Manual inspection of the list of upregulated genes suggests 
that multiple signalling pathways are affected upon FOX03 
activation (Figure ID). FOX03 induces an auto-regulatory loop 
as increased expression of PI3K pathway components is 
observed [IGF1R and IRS2). This likely results in increased 
PKB/AKT activation in the presence of growth factors 
ultimately inhibiting FOXO (Kloet and Burgering, 2011). In 
addition, RICTOR expression increases, which is an mTOR 
complex 2 component necessary for efficient FOXO phosphor- 
ylation and cytoplasmic retention (Brown et al, 2011). 
Importantly, activation of the PI3K/PKB axis downstream of 
FOXO activation has been suggested for both mammalian cell 
systems and model organisms (Ide et al, 2004; Puig and Tjian, 
2005; Hui et al, 2008), thus representing a conserved feature of 
FOXO biology. MAP kinase (MAPK) signalling is potentially 
inhibited upon FOX03 activation, in line with the observations 
in other studies (van den Heuvel et al, 2005) . MAPK activity 
negatively regulates FOX03 stability (Yang et al, 2008) and in 
this manner FOX03 could preserve its own expression. We also 
identified multiple ligands, receptors and other signalling 
components, involved in a variety of signal transduction 
pathways, including EGF, TGF, TNF, WNTand Rho signalling. 
This implies that FOX03 activation results in an overall change 
in the wiring or activity of several signal transduction 
networks. PMAIP1 (NOXA), a previously identified FOXO 
target involved in the regulation of apoptosis (Valis et al, 2011) , 
is also induced by FOX03 activation. Enrichment in cell-cycle- 
related genes can be observed for downregulated targets 
(Figure ID; Supplementary Table 1) as well as downregulation 
of genes with a role in mitochondrial biogenesis and function. 
This is in agreement with FOX03-induced inhibition of cell- 
cycle progression (Kops et al, 2002b) and mitochondrial 
function (Ferber et al, 2012). In general, the observed 
transcriptional changes suggest sensitization of the PI3K/ 
PKB axis, resulting in negative feedback on FOXO activation, 
inhibition of MAPK signalling and alteration of signalling 
networks, accompanied by inhibition of cell-cycle progression 
and mitochondrial function. 



FOX03 regulates target gene expression through 
RNAPII recruitment 

Transcription of eukaryotic genes is regulated at initial 
recruitment with additional regulation of pause release and 
elongation into the gene body (Nechaev and Adelman, 2011; 
Adelman and Lis, 2012). Using the RNAPII profiles, we sought 
to investigate whether FOX03-mediated target gene regulation 
occurs through regulation of RNAPII recruitment or pause 
release. If regulation of target gene expression would occur 
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Figure 1 Genome-wide FOX03-induced changes in RNAPII occupancy. RNAPII ChlP-seq was performed in DLD1-F3 cells containing 40HT-inducible F0X03 and 
DLD1 control cells. Profiles were generated in both cell lines without 40HT and 4 and 24 h upon addition of 40HT. (A) A representative example of RNAPII occupancy 
profiles. The genomic region surrounding the IRS2 gene locus is shown, for which RNAPII occupancy increases 2.8 (4 h) and 3.9-fold (24 h). V-axis values represent tag 
coverage per base per 10 6 sequenced reads. (B) Numbers of genes in RNAPII occupancy analysis. Total number of genes with detectable RNAPII occupancy and 
changed RNAPII occupancy specifically in DLD1-F3 cells, 4 and/or 24 h upon addition of 40HT. (C) Overlap of changes in RNAPII occupancy with changes in mRNA 
levels measured by micro-array gene expression analysis. Fold change in micro-array is distinguished in low (1 .3-2-fold) and high ( > 2-fold) induction. Hypergeometric 
test for overlap for the subset of RNAPII changed genes versus all RNAPII occupied genes (ns = P>0.05, *P<0.05, **P<5x 10~ 5 , ***P<5x 10~ 15 , 
****P< 5 x 1 0 ~ 30 ). (D) Overview of F0X03-regulated genes. Lead by processes identified in GO term analysis on changed genes, we manually categorized a subset of 
all genes changed in RNAPII occupancy upon FOX03 activation, confirmed by changes in mRNA levels. Guided by literature, we subcategorized genes involved in 
signal transduction, cell fate decisions and mitochondrial function (details in main text). 



through pause release, a shift in RNAPII occupancy from the 
TSS towards the gene body is expected. We performed meta- 
gene analysis of RNAPII occupancy over up- and down- 
regulated genes independently categorized by micro-array 
analysis. Upon induction, an overall increase in RNAPII 
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occupancy in DLD1-F3 profiles can be observed in the subset 
of upregulated genes (Figure 2A) . In addition, we observed no 
changes in the relative ratio between RNAPII signal around 
the TSS and in the gene body (travelling ratio, Figure 2B). 
For transcription factors involved in the regulation of pause 

Molecular Systems Biology 2013 3 



Genome-wide analysis of FOX03 gene regulation 
A Eijkelenboom et al 



DLD1 -40HT 
DLD1 4 h +40HT 
I DLD1 24 h +40HT 



IDLD1-F3 -40HT 
I DLD1-F3 4 h +40HT 
I DLD1-F3 24 h+40HT 



3 - 



O- 

< 




I I I 



Average upregulated gene 



TSS 



TES 



Gene 

(-/+1 kb) 5 ' UTR body 3 ' UTR (-/+1 kb) 
Location relative to gene annotation 



B 




10 



TSS TES 
DLD1-F3 -40HT 

• DLD1-F3 4h+40HT 

• DLD1-F3 24h+40HT 

Travelling ratio = 



3> < 



-10 



25 



50 



75 100 



% Of up-regulated genes 
(ranked) 



Figure 2 F0X03 induces target gene expression through recruitment of RNAPII. (A) Average RNAPII occupancy for upregulated genes independently classified by 
micro-array analysis. Genes were divided into 1 kb regions surrounding TSS and TES, the 5'UTR, 3'UTR and gene body. UTRs and gene body were normalized in size 
for every gene and average signal was determined for every position, excluding outliers (three highest and lowest values). Signal was normalized to 1 for the DLD1 -F3 
(-40HT) profile for all regions separately. (B) Travelling ratio for upregulated genes classified by micro-array analysis. Ratios for RNAPII signals (a= 0.3 kb up- and 
downstream of TSS, b= + 0.3 kb from TSS to annotated TES) are ranked by increasing ratio for DLD1-F3 untreated, 4 and 24 h treated separately. 



release, changes in this ratio have been observed upon 
alteration of transcription factor activity (Rahl et al, 2010). 
The preservation of overall distribution is consistent with 
FOX03 activation regulating the recruitment of RNAPII, 
without affecting pause release and the efficiency of early 
elongation. The increase around the TSS in 4 and 24 h time 
points is similar, suggesting that full induction of RNAPII 
recruitment is achieved within 4h of FOX03 activation. 
The downregulated subset does not show any prominent 
change in the average level of RNAPII (Supplementary 
Figure 2) . This could be due to the smaller overlap between 
micro-array and RNAPII experiments or different mechanisms 
of downregulation, which will be discussed later. 

Genome-wide binding profile of human FOX03 
reveals a transcriptional activating role for 
FOX03 binding 

To further understand the mechanism of FOX03 -induced 
transcription initiation, we analysed genome-wide binding of 
FOX03 by ChlP-seq. In total, we identified 9932 FOX03- 
binding sites (peak coordinates in Supplementary Table 2), 
similar to the amount of FOXOl-binding sites observed by 
others (Fan et al, 2010; Lin et al, 2010) . We confirmed 21 out of 
21 FOX03-binding sites with ChIP followed by qPCR 
(Supplementary Figure 3), showing reliability of the profile. 
We found increased FOX03 binding in several previously 
described target gene promoters analysed by ChlP-PCR and the 
majority by reporter assay studies (among others: MXI1 
(Delpuech et al, 2007), PIK3CA (Hui et al, 2008), PMAIP1 
(Valis et al, 2011), SOD2 (Kops et al, 2002a), PRDX3 (Chiribau 
et al, 2008), CCNG2 and KLF6 (Terragni et al, 2008)) although 
not all regions fulfil the peak calling criteria. 
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We performed de novo motif search to gain further insight 
into the sequence-requirements for FOX03 DNA binding and 
identified eight sequence motifs enriched in FOX03 bound 
regions compared with randomly selected regions. Identified 
motifs cluster close to the centres of FOX03-binding sites, are 
evolutionary conserved and enriched in both strong and also 
weaker FOX03-binding sites. Seven out of eight motifs were 
found to reliably associate with FOX03 peaks (Figure 3A; 
Supplementary Figure 4). Three of these motifs fulfil the 
Forkhead motif criteria (Furuyama et al, 2000; Biggs et al, 
2001), in line with direct DNA binding of FOX03 to these 
regions. Lower stringency in mapping increased the total 
fraction of Forkhead motif containing peaks to 61 % , showing a 
full or suboptimal Forkhead motif can be identified in the 
majority of peaks. In addition, another four identified motifs 
resemble known published motifs of AP-1 (Angel and Karin, 
1991), GATA (Ko and Engel, 1993), RUNX (Meyers et al, 1993) 
and SP-1 (Kadonaga et al, 1987), suggesting a possible role 
for these factors in FOX03 DNA-binding and target gene 
co-regulation. 

To further relate the location of peaks with respect to 
potential roles in gene regulation, we determined the distribu- 
tion according to the closest TSS for all FOX03 bound regions. 
The majority of identified binding sites is located outside of 
promoter regions in distal intergenic regions (> 5 kb from any 
annotated transcript), with a substantial fraction of binding 
sites found within gene bodies (33 % , Figure 3B) . Interestingly, 
the subset of distal FOX03-binding regions shows a high 
conservation around the peak centre (Figure 3C), indicating 
that these distal binding sites have biological relevance. 
To correlate changes in gene expression with FOX03 binding, 
we determined the distance from the TSS to the closest 
FOX03 peak for genes with changes in RNAPII occupancy 
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Figure 3 Genome-wide binding profile of human F0X03. (A) ChlP-seq analysis was performed to identify F0X03-binding regions in DLD1-F3 cells. De novo motif 
analysis revealed the presence of seven motifs enriched in the identified F0X03-binding regions. Relative enrichment to control regions and percentage of peaks 
containing the motif are shown. (B) Distribution of peak location with respect to annotated genes. Intragenic denotes within transcribed regions. In addition, the 
distribution relative to the closest TSS, compared with random regions, is shown. (C) Evolutionary conservation of FOX03 bound regions. Phastcons was used to 
produce base-by-base conservations scores of all FOX03 bound regions and distal intergenic bound regions located >5kb from the closest annotated transcript. 
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Figure 4 F0X03 binding correlates with upregulation of RNAPII occupancy. (A) Comparison of gene status in RNAPII occupancy and distance between TSS and 
closest FOX03 peak. The closest peak is determined by distance to the TSS only, irrespective of the location of other genes. The distribution of all RNAPII occupied 
genes is shown, as well as up- and downregulated (4 h, 24 h and both) genes (total number of genes in each category is indicated below x-axis). Hypergeometric test for 
number of genes with peak within 20 kb of subsets versus all (ns = P>0.05, *P<0.05, **P<5x 1(T 5 , ***P<5 x 1CT 15 , **"P<5 x 10 ). (B) Fold change in 
RNAPII occupancy (4 h) for groups of genes binned by increasing distance of gene TSS to closest F0X03 peak (group size is indicated above x-axis). The positive 
correlation is significant up to 20 kb (Mann-Whitney-Wilcoxon test, ns= P>0.05, *P<0.05, **P< 0.005, ***P< 0.0005, ""P< 0.00005). 



(Figure 4A). On average, upregulated genes are closer to 
FOX03 bound regions (up to 66% within 20 kb of the closest 
peak), compared with all genes (21%). Downregulated genes 
show a negative correlation (16%), but this is not significant 
due to the smaller differences and smaller number of 
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genes. Similar conclusions can be drawn from changes in 
mRNA levels (Supplementary Figure 5) . This suggests a pure 
activating role for FOX03 DNA binding in recruitment of 
RNAPII to target genes, which is in agreement with binding 
data of FOXO orthologues in C. elegans and Drosophila 
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(Schuster et al, 2010; Alic et al, 2011) and represents an 
evolutionary conserved feature of FOXO transcription factors. 



FOX03 binds to and activates enhancers 

We reasoned that FOX03 binding to distal regions might be 
relevant for transcriptional control, since we have observed 
binding of FOX03 to many distal regions with increased 
conservation. Genes with the closest FOX03 peak as distant as 
15-20kb still show a significant increase in RNAPII 
(Figure 4B), suggesting that FOX03 can activate target 
gene expression from greater distances. Enhancers are DNA 



elements located distant from promoter regions that can bind 
transcription factors, transcriptional co-activators and chro- 
matin regulators to affect transcription at promoters (Ong and 
Corces, 2011 ; Splinter and de Laat, 2011) . Previous studies have 
reported the presence of RNAPII and transcription at enhancer 
regions (De Santa et al, 2010; Kim et al, 2010; Ernst et al, 2011) . 
We therefore analysed RNAPII occupancy at distant FOX03 
bound regions and observed an increase in RNAPII levels upon 
FOX03 activation (Figure 5A) . 

To investigate the contribution of FOX03 -induced peaks to 
all intergenic RNAPII bound regions, we identified intergenic 
RNAPII peaks. A total of 2367 intergenic RNAPII peaks was 
found, similar to numbers from previous studies (De Santa 
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Figure 5 F0X03 recruits RNAPII to intergenic regions. (A) RNAPII levels at F0X03 bound intergenic regions. Intergenic regions are located over 5 kb from any known 
transcript or known H3K4me3 peak. RNAPII signal is specifically increased in DLD1-F3 4h and 24 h following 40HT addition, with no change in DLD1 control cells 
(median within box, P-values from Mann-Whitney-Wilcoxon test, treated versus untreated). (B) RNAPII levels only increase in FOX03 bound intergenic RNAPII peaks. 
Overlap between intergenic RNAPII peaks and intergenic FOX03 peaks or matched random regions is determined. Overlap is higher as expected by chance with 
FOX03 peaks compared with matched random regions (Pearson's z -test with Yates' continuity correction). Boxplot shows levels of RNAPI I occupancy only increases in 
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the FOX03 peak centre. Regions are sorted by descending RNAPII signal in DLD1-F3, 4h. Signal intensities are normalized for sequencing depth. Maximum colour 
intensity corresponds to the 0.98th percentile of the normalized read density distribution. Graphs represent average RNAPII signal of the same 8 kb regions, line colours 
are indicated below heat maps. 



6 Molecular Systems Biology 2013 



© 2013 EMB0 and Macmillan Publishers Limited 



Genome-wide analysis of FOX03 gene regulation 
A Eijkelenboom ef a/ 



et al, 2010). 637 peaks overlap with FOX03 peaks (Figure 5B) 
accompanied by a specific increase in RNAPII occupancy in 
FOX03 bound RNAPII peaks only, thereby excluding global 
changes in RNAPII ChIP efficiency upon FOX03 activation. In 
an alternative approach, we created three classes of intergenic 
RNAPII peaks after FOX03 induction based on occupancy 
status: increased, decreased and unchanged (Supplementary 
Figure 6A and B). The increased class shows the greatest 
overlap with FOX03 peaks and, independently of the FOX03 
ChlP-seq, we identified the Forkhead motif by de novo motif 
search in these regions. 

To further analyse the RNAPII binding profiles and changes, 
we generated heat maps covering 8 kb surrounding all FOX03 
bound distal intergenic regions and plots of the average signal 
intensity at these locations (Figure 5C). Prior to FOX03 
binding, the RNAPII intensity is high surrounding the FOX03 
peak centre, showing that RNAPII is already present in the 
majority of FOX03 bound regions preceding induction. Upon 
FOX03 activation, the largest RNAPII occupancy change can 
be observed directly at the FOX03 peak centre and spreading 
up to several kb surrounding the FOX03-binding site. 

Combinations of modifications of the histone tails reflect the 
function and status of a locus (Zhou et al, 2011) . The canonical 
chromatin signature of distant enhancers is the presence of 
Histone H3 Lysine 4 monomethylation (H3K4mel), while 
trimethylation (H3K4me3) is generally absent (Heintzman 
et al, 2007, 2009). Enhancers are also marked by the presence 
of H3K27 acetylation (Ac), associated with active enhancers 
(Creyghton et al, 2010). We performed ChlP-qPCR with 
antibodies recognizing specific histone modifications and 
show five FOX03 bound regions are high in H3K4mel, while 
H3K4me3 is low, with an increase in H3K27AC following 
FOX03 activation (Figure 6A) . We probed regions surrounding 
the FOX03 peaks and could confirm the local enrichment 
of H3K4mel and local increase of H3K27Ac (Supplementary 
Figure 6C). The presence of H3K27Ac prior to FOX03 acti- 
vation varies between tested regions. Similar variation can be 
observed for RNAPII occupancy (Figure 5C; Supplementary 
Figure 6D) , likely reflecting differences in activity. In addition, 
four out of five sequences were responsive to FOX03 induction 
in a luciferase-reporter assay (Figure 6B) . This activation was 
orientation independent, a hallmark of enhancers (Banerji 
et al, 1981). Together, we show that distant FOX03-binding 
sites are marked with a combination of histone modifications 
specific for enhancers and can obtain a more active signature 
following FOX03 activation with increased histone acetylation 
and RNAPII occupancy. 



FOX03A3-ER faithfully mimics endogenous 
FOX03 

We made use of an inducible system to study direct effects 
of FOX03 activation. Although considered to act neutral, we 
wished to exclude any influence from the ER-tag. Therefore, 
we generated DLD1-DBD cells containing a fusion of the ER 
moiety to only the FOX03 DNA-binding domain (DBD) and 
Nuclear Localization Signal (NLS) (Supplementary Figure 7A) . 
Upon activation with 40HT, this fusion protein should bind to 
the same regions as full-length FOX03, but as it lacks FOX03 
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Figure 6 FOX03 binds to and activates enhancers. (A) ChlPs on 
DLD1-F3 cells (not treated and 2h following 40HT treatment) with antibodies 
recognizing FOX03, H3K4me1, H3K4me3 and H3K27Ac. Immunoprecipitated 
DNA was analysed by qPCR for five FOX03 bound distal regions (closest 
gene-distance-location relative to TSS: A= DL/SP5-25kb-downstream, 
6= 7GFBP2-128kb-upstream, C= PDGFRA-5 kb-downstream (intron), D= 
HMGA129kb-downstream (intron), E= KC6-210kb-upstream) and three 
unbound regions (GD = GAPDH promoter, MB = Myoglobin locus and NC3 = 
negative control). GD is within the GAPDH promoter region as a positive control 
for histone ChlPs. Percentage of input was calculated and enrichment relative to 
MB control region is shown for three to five independent biological replicates 
(mean with s.e.m., P-value from paired f-test, treated versus untreated, 
*P<0.06, **P<0.05, ***P< 0.005). (B) Luciferase assay on constructs 
containing a minimal promoter with upstream insertion of the same five FOX03 
bound regions. Orientation specificity was tested by analysing both orientations 
(F and R) of the same sequence. TK-renilla was co-transfected as internal control 
and luciferase/renilla values were determined. Fold change after 8h 40HT 
versus untreated were calculated for DLD1 and DLD1-F3 cells (mean of three 
independent biological replicates, s.e.m., P-value from paired f-test, DLD1-F3 
versus DLD1 , *P< 0.05, **P< 0.005, ***P< 0.0005). Source data for this figure 
is available on the online supplementary information page. 
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transactivation domains, this will allow us to observe any side 
effects induced by the ER-tag. We first confirmed binding of the 
DBD fusion to six regions previously identified to bind full- 
length FOX03A3-ER (Supplementary Figure 7B). The binding 
of the ER-DBD did not induce H3K27 acetylation at these 
regions (Supplementary Figure 7C) . To observe effects of the 
DBD fusion on RNAPII occupancy, we generated genome- wide 
profiles of ER-DBD binding and RNAPII changes upon 40HT 
treatment in DLD1-DBD cells. We found substantial overlap 
between ChIP profiles (e.g., 18% of DBD peaks overlaps with 
called FOX03A3-ER peaks) and in the vast majority of not 
overlapping peaks we found a clear signal which was below 
peak detection limit (Supplementary Figures 8 and 9A). 
The DBD fusion protein binds to intergenic FOX03 bound 
regions, but we did not observe any changes in RNAPII 
occupancy at these regions or at intergenic DBD bound regions 
(Supplementary Figure 9) . 40HT treatment in DLD1-DBD cells 
also does not induce any similar changes as FOX03A3-ER in 
gene expression, judged by both RNAPII occupancy 
(Supplementary Figure 10) and mRNA levels (Suppleme- 
ntary Figure 11). From this we conclude that the ER-tag acts 
neutral and that the FOX03A3-ER fusion protein requires the 
FOX03 transactivation domains to mediate the observed 
effects on histone acetylation and RNAPII occupancy. 

We employed PI3K/PKB inhibition to induce and study 
endogenous FOX03 activation. First, we confirmed binding of 
endogenous FOX03 to the same six regions described above 
and PKB inhibition resulted in increased H3K27Ac at these 
locations, resembling the increase observed upon activation of 
exogenous FOX03 (Supplementary Figure 7) . Second, we also 
generated genome-wide binding profiles of endogenous 
FOX03, in which we found considerable overlap between 
endogenous and exogenous ChIP profiles (Supplementary 
Figure 8), revealed by overlap in called peaks (34% of 
endogenous FOX03 peaks overlaps with called FOX03A3-ER 
peaks) and increased signal around peak centres from both 
data sets (Supplementary Figures 8B and 9A) . In addition, to 
compare the above-described effects of exogenous FOX03 
activity with the effects of PKB inhibition, we generated 
RNAPII profiles in DLD1-F3 cells comparing addition of 40HT 
and PKB inhibitor. At the distal intergenic endogenous 
and exogenous FOX03 bound locations identified above, 
PKB inhibition mimics the increase in RNAPII occupancy 
upon activation by 40HT (Supplementary Figure 9) . Next, we 
calculated RNAPII occupancy changes at gene coding regions 
and compared these with above-described results (Figure 1). 
We identified the majority of genes (Supplementary Figure 10) 
to again be changed by 40HT treatment: i.e., over 51% (0.5- 
fold cutoff) or 83% (0.2-fold cutoff) of genes previously 
identified to be upregulated by 4- and 24-h treatment with 
40HTwere again identified to be upregulated. In addition, PKB 
inhibition induced similar changes with less fold induction, 
showing 17% (0.5-fold) and 60% (0.2-fold) of these same 
genes to be upregulated. Finally, we confirmed induction in 
mRNA levels for six genes upon activation of exogenous 
FOX03 and PI3K/PKB inhibition (Supplementary Figure 11), 
repeatedly showing clear induction but slightly lower fold 
change compared with FOX03A3-ER. 

Our analysis shows that upon PI3K/PKB inhibition, endo- 
genous FOX03 binds to the same regions as FOX03A3-ER. We 



also show concomitant similar changes on RNAPII profiles and 
histone acetylation. Generally, changes induced by PKB 
inhibition are somewhat smaller relative to activation of 
FOX03 A3-ER. Likely the FOX03A3-ER system provides a more 
robust profile as FOX03 activation through PKB inhibition still 
requires dephosphorylation of FOX03, which likely results 
in a slower response. Alternatively, PKB inhibition might not 
induce the same level of constitutive activation like the A3 
mutant does, since PKB inhibition will likely also affect 
signalling to other transcription factors that for example may 
enhance or reduce FOXO only effects. Thus, isolated activation 
of FOX03 in our inducible system could also lack activation of 
other factors possibly cooperating with FOX03 in target gene 
regulation. The great overlap in 40HT and PKB inhibition 
induced changes (Supplementary Figure 10C), however, 
suggests that generally FOX03 requires cooperation with 
activators already active under conditions studied. 



FOX03 bound enhancers communicate with target 
gene promoters through pre-existing DNA loops 

We initially employed the ER inducible system because we 
aimed to identify direct effects of FOX03 activation. Unlike 
inhibition of the PI3K/PKB pathway, this system allows us to 
relate all observed changes specifically to activation of only 
FOX03. We therefore sought to investigate a possible relation 
between the increase of RNAPII at intergenic FOX03 bound 
regions and FOX03-induced target gene regulation. Transcrip- 
tion of RNAPII occupied enhancer regions has been reported to 
precede the activation of adjacent protein-coding genes 
(De Santa et al, 2010). Indeed, we could observe a positive 
correlation between the RNAPII increase at FOX03 bound 
regions and upregulation of the adjacent genes (Figure 7A; 
Supplementary Figure 12A), which could be explained by 
communication between FOX03 bound regions and the 
neighbouring target gene. To examine whether FOX03 bound 
regions can physically loop towards the TSS of a target gene, 
we performed high-resolution chromosome conformation 
capture analysis followed by sequencing (4C-seq) (Splinter 
et al, 2011). 4C can be used to generate a map of the spatial 
interactions of a selected locus. Based on visual inspection of 
regions around FOX03-induced target genes, we selected 
candidate regions with increased RNAPII occupancy at both 
FOX03 bound putative enhancer regions and target genes. 
Primers were designed for seven view points, consisting of 
four intergenic FOX03 bound regions and three TSSs from 
induced genes (Supplementary Table 3) . Figure 7B shows the 
region around the KLF5 locus, for which the RNAPII 
occupancy increases (1.9-fold after 4h). 4C primers were 
designed around a FOX03 bound region 18kb upstream from 
the TSS for which the RNAPII occupancy also increases 
(2.1 -fold after 4h). Typical for a 4C profile, the frequency of 
captured genomic regions decreases as a function of the 
distance to the viewpoint. However, capture frequency greatly 
increases around the TSS of KLF5, indicating looping of 
the FOX03 bound intergenic region towards the TSS region. 
The signal also increases slightly around other FOX03 bound 
regions, hinting at the close proximity of multiple enhancers to 
the target gene TSS. The 4C profile is highly similar before and 
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after FOX03 induction, suggesting that FOX03 binds to 
pre-existing loops and does not induce reorganization of the 
spatial environment of the locus. For three other regions, we 
find evidence of pre-existing loops between FOX03 distant 
binding sites and target genes (Supplementary Figure 12). 
For the IRS2 gene, where we look from the TSS, we see looping 
to two upstream FOX03-binding sites. It is important to note 
that the closest (and strongest) FOX03 peak is not engaged in 
an interaction with the TSS of IRS2. Together, these results 
suggest pre-existing communication between several FOX03 
bound distal regions and regulated target genes. 



Discussion 

Using genome-wide approaches we have analysed how FOX03 
regulates target gene expression. To this end, we measured 
FOX03 binding to genomic DNA and changes in RNAPII 



occupancy following FOX03 activation. This analysis reveals 
several features of FOX03-mediated gene regulation. 

First, binding of FOX03 to gene regulatory regions activates 
transcription and gene expression, whereas repression of gene 
expression following FOX03 binding generally occurs through 
DNA-binding independent or indirect mechanisms. FOX03- 
induced MYC-mediated gene repression possibly entails an 
important mechanism by which FOXO indirectly represses 
gene expression. We observe an increase in MXI1 expression, a 
previously reported FOX03-induced negative regulator of MYC 
function (Delpuech et al, 2007). The MYC DNA-binding 
sequence is enriched in promoter regions of downregulated 
genes (CACGTG, (Blackwell et al, 1990), Supplementary 
Figure 13). This is in agreement with other studies that show 
an intricate and extensive inverse correlation between FOXO 
and MYC function (Kloet and Burgering, 2011). In addition, 
induction of cell-cycle inhibitor CCNG2 could be responsible 
for downregulation of cell-cycle genes. Other cell-cycle 
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regulators previously reported to be regulated by FOXO 
(van der Vos and Coffer, 2011) were not changed above the 
thresholds used in this study, but possibly still contribute. 
Alternatively, the smaller overlap between RNAPII and mRNA 
levels in downregulated genes could reflect changes in mRNA 
stability rather than transcriptional changes. Although the 
activating role for FOX03 DNA binding in gene expression is 
conserved and reported for C. elegans and Drosophila FOXO 
orthologues as well (Schuster et al, 2010; Alic et al, 2011) , some 
examples of FOXO-mediated suppression of target gene 
expression have been reported. This includes repression of 
CCND1 (Schmidt et al, 2002), ID1 (Birkenkamp et al, 2007) 
and, during revision of this manuscript, FOX03-mediated 
suppression of several antiviral genes in macrophages was 
reported (Litvak et al, 2012). Although our data set clearly 
shows FOX03 DNA binding is correlated with transcriptional 
activation only, we do not wish to exclude the possibility that 
FOX03 activation induces repression of individual genes. 

Second, the changes observed in the RNAPII profile upon 
FOX03 activation are consistent with FOX03 regulating 
transcription initiation rather then pause release and elonga- 
tion. This is in contrast to yeast forkheads, which have been 
reported to regulate RNAPII elongation (Morillon et al, 2003). 
We observe FOX03 binding within the coding region of genes 
for roughly one-third of the peaks. In general, we can detect a 
slightly higher signal in the FOX03 ChlP-seq throughout 
transcribed genes, but it is unclear whether this reflects 
genuine FOX03 binding or a nonspecific increase in back- 
ground signal. Binding of yeast forkheads and Drosophila 
FOXO within transcribed regions has been observed as well 
and has been proposed to be necessary for elongation 
(Morillon et al, 2003; Alic et al, 2011). Since our findings 
argue against a role for FOX03 in pause release and elongation, 
the relevance of FOX03 binding within genes is unclear. 

Generally, pausing of RNAPII after initial recruitment seems 
specifically enriched at genes in signal-responsive pathways 
(Adelman and Lis, 2012) possibly allowing a quick response in 
gene expression upon stimulation. In light of this, we consider 
it surprising to only find evidence for the regulation of 
initiation by FOX03 activation. 

Third, the FOX03 ChlP-seq resulted in the identification of 
nearly ten thousand FOX03-binding sites. In addition to the 
Forkhead motif, we identified four motifs enriched in 
FOX03 bound regions, resembling known published motifs 
for AP-1, GATA, RUNX and SP1. ChlPs on FOXO orthologues 
in Drosophila adults and C. elegans have identified AP-1 
(C. elegans only) and GATA (both) binding sequences in a 
subset of FOXO bound regions as well, suggesting that the co- 
occurance is evolutionary conserved (Schuster etal, 2010; Alic 
et al, 2011). We could, however, not find any correlation 
between motif presence and RNAPII changes. From our data 
set it is therefore unclear to what extent these motifs, and the 
proteins that bind these sequences, influence FOX03-induced 
target gene expression. 

Pioneer transcription factors are transcription factors 
thought to be capable of initial binding to regulatory 
sequences, allowing binding of other factors by opening 
compacted chromatin, initiating the events that will ultimately 
result in transcriptional activation. Forkhead transcription 
factor FOXA1 has been studied extensively for its pioneering 
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activity (reviewed in Zaret and Carroll, 2011). FOXOl has been 
shown to similarly open compacted chromatin in vitro (Hatta 
and Cirillo, 2007; Hatta et al, 2009) suggesting that it may 
function as a pioneer factor as well. Reports on correlation 
of changes in histone modifications and FOXA1 activity are 
inconsistent, either reporting the presence of H3K4mel 
upstream of, and importance for, FOXA1 binding and/or 
recruitment (Lupien et al, 2008) or induction of methylation 
after FOXA1 binding to regions lacking H3K4mel (Serandour 
et al, 2011) . Based on this it is unclear whether the higher level 
of H3K4mel at several FOX03-binding sites conflicts with the 
hypothesis of FOX03 as a pioneer. The genome-wide RNAPII 
occupancy profiles revealed the local presence of RNAPII prior 
to FOX03 binding in the majority of studied regions, 
suggesting that in general FOX03-binding regions are already 
in a more open and potentially active state. This would argue 
against FOX03 acting as a pioneering factor in these regions 
under the conditions used in this study, without excluding that 
FOX03 acts as a pioneer in a subset of bound regions that 
lack RNAPII signal before FOX03 binding. GATA has been 
described to function as a pioneering factor as well (reviewed 
in Zaret and Carroll, 2011), and RUNX1 was shown to be 
transiently required for chromatin unfolding for initial 
transcription (Hoogenkamp et al, 2009). The enrichment of 
GATA and RUNX binding motifs within FOX03 bound regions 
might therefore indicate these to act as pioneering factors that 
thereby cooperate with and/or precede FOX03 binding. 
Alternatively, FOXA1 may act as a pioneering factor for FOX03 
since FOXA1 binds the same Forkhead motif as FOX03. 
However, for the same reason this possibility can therefore 
also not be addressed by our data through analysis of motif 
presence and would require detailed time course analysis of 
genome binding for multiple Forkhead members, including 
FOXA1 and FOX03. 

Fourth, we found a function for distal FOX03 bound regions 
in target gene regulation. These regions are marked by the 
presence of canonical enhancer specific histone modifications 
and display a more active signature upon FOX03 activation, 
including an increase in H3K27Ac and recruitment of RNAPII. 
Direct binding between p300 histone acetyltransferase and 
FOXOs has been observed in many studies (Calnan and 
Brunet, 2008). Since p300 has been linked to transcription 
initiation (Kraus et al, 1999) and binding to enhancers has 
been reported (Heintzman et al, 2007; Xi et al, 2007; Visel et al, 
2009), recruitment of p300 by FOX03 could be responsible for 
the observed increase in H3K27Ac at enhancers and RNAPII 
recruitment to enhancers and target genes. 

Genome-wide studies for FOXO binding in Drosophila and C. 
elegans have not reported binding to distant sites, perhaps in 
part because these genomes are much smaller. Potentially, 
FOXO regulation of transcription from distant sites might be an 
acquired feature of mammalian FOXO biology not present in 
insects and nematodes. Other mammalian studies have 
identified FOXO binding in promoter regions and also at more 
distal locations, possibly at enhancers (Fan et al, 2010; Lin etal, 
2010; Chen et al, 2011), but to what extent these contribute to 
regulation of gene expression was unclear. To our knowledge, 
no studies have reported looping of FOXO occupied regions to 
promoters of target genes. Over the past years, several proteins 
including GATA and CTCF, have been reported to be involved 
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in looping in transcriptional regulation (Sexton et al, 2009). 
This could be an alternative explanation for the identification 
of the GATA motif in FOX03 bound regions. In addition, the 
CTCF motif was found enriched in Drosophila FOXO bound 
regions (Alic et al, 2011) . However, regulation and formation of 
enhancer-promoter looping is complex and establishing any 
involvement of these factors in FOXO target gene regulation 
requires thorough investigation. We observed that several 
FOX03 bound regions might fold to one promoter. Regulation 
of target genes through several regulatory regions simulta- 
neously could be necessary for tightly controlled gene 
expression, for example, through setting a threshold. This 
might therefore also explain why the number of FOX03 bound 
regions greatly outnumbers the sum of regulated genes. 
Alternatively, many FOX03 bound regions might not be 
functional in gene regulation in a given tissue or in general. 
In C. elegans a similar observation of DAF-16-binding sites 
outnumbering DAF-16-regulated genes has led to the sugges- 
tion that many sites are 'parking sites' occupied to ensure rapid 
DAF-16 transcriptional response when required (Schuster 
et al, 2010). Although we are not aware of reports showing 
chromatin looping being present and functionally relevant in 
C. elegans, our observation that several FOX03-binding sites 
may co-localize in three dimensions clearly indicates an 
alternative possibility besides simple 'parking'. 

The pre-existing nature of the enhancer-promoter loops and 
the cell type specificity of enhancers and chromatin architec- 
ture (Ernst et al, 2011; Sanyal et al, 2012) suggest that cells are 
poised for specific behaviour. Variation in overall chromatin 
architecture between different systems may therefore at least 
in part explain the previously mentioned context-dependent 
nature of FOXO activation. Assessing chromatin architecture 
around differently responding target genes may prove 
necessary to determine the regulatory regions responsible for 
FOXO-induced changes in gene expression. 

Our analysis has provided insight into the molecular 
mechanisms by which FOX03 induces target gene 
regulation. We have revealed a function for FOX03 outside 
of promoter regions. We have observed binding of FOX03 to 
enhancers accompanied by the acquirement of a more active 
signature. Finally, our data suggest importance of pre- 
existing enhancers and chromatin architecture for FOX03 
activated transcription and context-dependent gene regulation 
(Figure 7C). 



Materials and methods 
Cell culture 

DLDl, DLD1-F3 and DLD1-DBD cells were grown on RPMI-1640 
medium with 10% FCS and standard supplements. 4-Hydroxy- 
tamoxifen (40HT) from Sigma was dissolved in ethanol and added 
to cells at a final concentration of 1 \iM. Cells were treated with PI3K 
inhibitor LY 294002 (Enzo Life Sciences) at a concentration of 10 uM 
and PKB inhibitor VIII (Santa Cruz) at a concentration of 10 uM. 
pcDNA3-ER-DBD-NLS was generated by amplification of DBD-NLS 
domains of FOX03, which were amplified from FOX03A3 plasmid 
using forward primer GCGCGCGGATCCGCTGGGGGCTCCGGGCAG 
and reverse primer GCGCGCCTCGAGTCAGGCTGCGCGGCCACGGCT, 
digested with BamHl and Xhol and ligated into pcDNA3-ER vector. 
DLD1-DBD cells were generated by transfection with fugene HD 
(Promega), selected with G418 (Santa Cruz, 500ng/ml) and clonal 



lines were generated by limited dilution. Protein levels were 
determined by SDS-PAGE and western blot, with anti-ER (Santa Cruz, 
MC-20) and anti-tubulin (Calbiochem, CP06). 

RNAPII chromatin immunoprecipitation and 
sequencing 

Approximately 20 x 10 s DLDl and DLD1-F3 cells grown in the absence 
or presence of 40HT for 4 and 24 h were used for ChlP-seq procedure. 
ChIP was performed according to previously described procedure 
(Hatzis et al, 2008; Mokry et al, 2010), details in Supplementary 
methods, using anti-RBPl (PB-7C2) antibody (Euromedex). Chroma- 
tin was additionally sheared, end-repaired, sequencing adaptors were 
ligated and the library was amplified by LMPCR. After LMPCR, the 
library was purified and checked for the proper size range and for 
the absence of adaptor dimers on a 2% agarose gel and sequenced on 
SOLiD/AB sequencer to produce 50-bp long reads. Sequencing reads 
were mapped against the reference genome (hgl9 assembly, NCBI 
build 37) using the BWA package (Li and Durbin, 2009). Nonuniquely 
placed reads were discarded. If more than five reads mapped to the 
same location and strand, only five reads were used for further 
analysis. For duplicate experiments, two times 20 x 10 s DLD1-F3 cells 
were grown in the absence or presence of 40HT or PKB inhibitor 
for 4h. For profiling in DLD1-DBD cells, 20 x 10 s cells were grown in 
the absence or presence of 40HT for 4 h. Sequencing reads were again 
mapped against the reference genome and nonuniquely placed reads 
were discarded. 



Quantification of RNAPII occupancy 

To set gene expression from RNAPII ChlP-seq data, we counted the 
number of the sequencing tags aligned to annotated transcript 
coordinates (Mokry et al, 2012). To avoid transcripts with zero 
mapped tags to interfere with logarithmic transformation of read 
counts, one read per every 10 million sequencing tags was added to 
each transcript. Raw read counts were normalized to the transcript 
length (from TSS to TES) and sequencing depth and quantile 
normalized. All samples are presented as quantile normalized read 
counts per transcript per kb of transcript per million sequencing tags 
(NRPKM), log2 scale. Expressed gene is gene with an expression 
higher than -1.5 NRPKM in at least one condition within an 
experiment. Differentially transcribed genes were set as genes with 
at least 0 . 5 (log2) fold NRPKM change in DLD 1 -F3 4 and 24 h compared 
with untreated, while fold change in DLDl cells is less than 0.2 (log2). 
Changes in RNAPII occupancy in DLDl samples were used to estimate 
noise and determine cutoff values (Supplementary Figure 1). 
If differences in fold change between annotated transcripts of the 
same gene occur, the highest fold change was used for further analysis. 
To identify differentially transcribed genes in replicate experiments 
using 4 h of 40HTor PKB inhibition, average NRPKM was determined 
prior to the calculation of fold changes. For comparison between 
data sets (Supplementary Figure 10), all genes below - 1.5 NRPKM 
expression threshold in either one of the experiments or changed over 
0.2-fold in DLDl cells following 40HT treatment were excluded from 
analysis. 



GO term analysis 

GO term analysis was performed using Database for Annotation, 
Visualization and Integrated Discovery (DAVID) . Regulated genes were 
determined by a 0.5 -fold (log2) change in RNAPII occupancy at 4 
and/or 24 h in DLD1-F3 only, for which changes in mRNA levels are 
confirmed in micro-array gene expression analysis. All RNAPII 
occupied genes were used as a background list. 

Determination of intergenic RNAP II peaks 

Cisgenome software package (li et al, 2008) was used for the 
identification of RNAPII peaks from the ChlP-seq data. Regions within 
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5kb of a known transcript (RefSeq) or H3K4me3 peak from a 
published data set in DLD1 cells [Tanimoto et al, 2010) are excluded. 
Regions larger than 3 kb in size were discarded as putative 
not-annotated genes. Only peaks identified in at least two out of four 
samples (DLD1 Oh, DLD1-F3 Oh, 4h and 24 h) were included. 



Online access data sets 

ChlP-seq and 4C raw data from this publication have been deposited 
at GEO with the accession numbers GSE35486 (ChlP-seq) and 
GSE36835 (4C). 



Chromatin immunoprecipitation of FOX03 and 
histone modifications 

DLDl, DLD1-F3 and DLD1-DBD cells were grown in the absence or 
presence of 40HTor PKB inhibitor. For ER and FOX03 ChlP-seq, 40 x 10 6 
DLDl and DLD1-F3 cells were used, while 10 x 10 s DLD1-F3 cells were 
used for histone modifications. IPs were performed for 16 h with the 
following antibodies: 1 ug of H3K4mel, H3K4me3 and H3K27AC and 
total H3 from Abeam (ab8895, ab8580 and ab4729, abl791), 5 ug of oeER 
or aFOX03 (Santa Cruz, MC-20 and H144) and 5 ug of normal rabbit IgG 
(Santa Cruz) . ChlPs were performed according to a previously described 
protocol from the Brunet lab (Renault et al, 2011), details in 
Supplementary methods. Sequencing library preparation, sequencing 
and mapping were performed for FOX03 samples as for the RNAPII 
ChlP-seq. Multiple reads mapping to same location and strand have been 
collapsed to a single read and only uniquely placed reads were used for 
peak-calling. Cisgenome software package (Ji et al, 2008) was used for 
the identification of binding peaks from the ChlP-seq data and further 
analysis. Cisgenome 1.1 was used for calling initial DLD1-F3 aER data 
set: -c0.2, -wlOO, - s 25, - p 0.5, - g 200 -1 225, DLDl aER sample 
was used as background, For all other data sets Cisgenome 2 with 
settings: — e 50, — maxgap 200, minlen 200; was used; for endogenous 
FOX03 ChlPs IgG was used as background, for RNAPII human input 
DNA from (Mokry et al, 2010). Conservation was determined by 
Preformatted PhastCon scores downloaded from Cisgenome. qPCR 
analysis was performed with IQ SYBR-Green mix (Bio-Rad) , see primer 
sequences in Supplementary Table 3 . 



De novo motif discovery 

Motif analysis was done using CisModule function (Zhou and Wong, 
2004) incorporated in Cisgenome sofware, using the following 
parameters: motif number K= 15, mean motif length Lambda = 10, 
maximal motif length allowed =18, initial motif length =10, initial 
module size D = 3.0, module length =100, order of background 
Markov chain = 3, MCM iteration = 500. For lower stringency in 
Forkhead motif mapping, r=50 instead of r=500 was used ( — r is 
defined in the Cisgenome utility: motifmap_matrixscan_genome). 
A combination of Cisgenome functions, custom PERL and R scripts 
was used for additional data analysis. Similarity search with known 
transcription factor binding motifs was performed using TOMTOM 
motif comparison tool (Gupta et al, 2007) . 



Luciferase assay 

DLDl and DLD1-F3 cells were transfected using Fugene HD (Promega), 
with TK-renilla and pGL3-promoter vector without insert or with inserts 
generated by PCR and digestion with Kpnl (primers in Supplementary 
Table 3) from genomic DLD1-F3 DNA. Lysis and subsequent determina- 
tion of luciferase and renilla counts was performed 48 h after transfection 
using the Dual Luciferase Reporter Assay (Promega) . 



4C 

4C analysis was performed with 20 x 10 6 DLD1-F3 cells cultured in the 
absence or presence of 40HT for 4 h as described previously (Splinter 
et al, 2011), but with Dpnll and Csp digestion. PCR was performed on 
10 times 100 ng of 4C template for each primer pair and pooled prior to 
purification. Primers are listed in Supplementary Table 3. 

Statistical analysis 

Analysis were performed in R, Excel or Graph Pad. The details of tests 
used are given in figure captions. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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